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RESUME 


Ce  rapport  presente  une  generalisation  de 
l'estimateur  de  luenberger  pour  les  systemes  stochas- 
tiques  et  deterministes  discrets  et  a delai.  Dans  de 
tels  systemes,  il  est  possible  que  les  signaux  detec- 
tes  ne  contiennent  que  du  bruit  et  que  seule  la  pro- 
bability que  ces  cas  surviennent  est  connue  de  1'es- 
timateur.  Cette  generalisation  contient  les  cas  par- 
ticulars du  filtre  optimal  de  Kalman,  du  filtre  de 
Nahi  et  de  l'estimateur  d'ordre  reduit  de  Tse  et 
Athans,  et  Leondes  et  Novak.  L'utilisation  de 
cette  generalisation  est  illustree  par  un  exemple 
de  poursuite  radar.  (NC) 


ABSTRACT 


This  report  presents  the  generalization  of 
Luenberger's  observer  to  discrete  time-delayed  deter- 
ministic and  stochastic  systems,  when  the  detected 
signals  may  contain  noise  alone,  and  only  the  proba- 
bility of  occurrence  of  such  cases  is  known  to  the 
estimator.  The  solution  encompasses  the  special  cases 
of  Kalman's  optimal  filter,  Nahi's  filter,  Tse  and 
Athans  and  Leondes  and  Novak's  reduced -order  obser- 
vers. The  design  of  a reduced-order  observer  is 
illustrated  using  a radar-tracking  problem.  (U) 
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1.0  INTRODUCTION 

Many  practical  problems  in  navigation,  weapon  delivery  and 
tracking  systems  require  the  use  of  the  Kalman  filtering  technique  to 
minimize  the  estimation  errors  of  certain  physical  variables.  In 
practical  realization,  however,  the  complexity  of  the  optimal  filtering 
technique  may  unduly  burden  the  computer  memory  and  computational  capa- 
bility. To  circumvent  this  problem,  there  has  been  much  interest  in 
the  construction  of  optimal  and  suboptimal  reduced-order  observers, 
along  the  line  of  Luenberger's  observer  [l].  One  group  of  papers 
[2-3]  essentially  obtained  the  conditional  mean  estimate  of  the  state, 
while  minimizing  the  order  of  the  filter.  The  other  group  [4-6]  con- 
structed the  observer  estimator  by  constraining  the  memory  size  of  the 
. observer.  And  they  all  made  the  classical  assumption  that  the  measure- 

ment signals  received  always  contained  the  state  variables. 

In  many  realistic  military  applications,  the  presence  of 
the  state  vector  may  be  associated  with  a non-zero  probability  of  false 
alarm  due  to  fading  or  enemy  interference.  In  the  case  that  all  the 
measurement  outputs  are  corrupted  by  additive  white  noise,  one  can  use 
the  filtering  algorithm  of  Nahi  [7].  But  to  the  author's  knowledge, 
optimal  and  suboptimal  observer-estimators  have  not  been  derived  in  the 
literature. 

Therefore,  this  report  is  concerned  with  the  generalization 
and  extension  of  Luenberger's  observer  to  discrete  time-delayed  deter- 
ministic and  stochastic  systems,  when  the  signal  detection  process  is 
associated  with  a non-zero  probability  of  false  alarm.  The  derivation 
makes  use  of  the  matrix  minimum  principle  [10]  to  minimize  the  error- 
variance  estimation  criterion. 

Section  2 presents  an  observer;  its  structural  design  can 

I’  considered  as  the  extension  of  Aoki  and  Huddle's  [4]  fundamental 

structure  in  the  domain  of  linear  deterministic  discrete  time-delayed 
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systems  under  uncertain  observations.  Section  3 applies  the  matrix 
minimum  principle  to  yield  the  minimum  variance  estimate  of  the  state. 
This  section  considers  the  situation  where  some  of  the  components  of  the 
measurement  vector  are  noise-free,  and  the  development  is  parallel  to 
that  of  Yoshikawa  and  Kabayaski  [8].  Section  4 presents  a class  of 
general  reduced-order  observers  under  uncertain  observation,  based 
upon  the  experience  of  Leondes  and  Novak  [5,6].  The  results  can  be 
shown  to  encompass  the  special  cases  of  Kalman  [9],  Nahi  [7]  and 
Leondes  and  Novak  [5,6].  A practical  radar-tracking  problem  is  pre- 
sented for  the  purpose  of  illustration. 

The  work  was  performed  at  DREV  duTing  the  fiscal  year  1977/78 
under  PCN  21J05,  "Control  Concepts". 


2.0  DETERMINISTIC  OBSERVER 


Consider  a general  class  of  linear  time-varying  discrete 

systems  involving  multiple  delays: 

N M 

x(k+l)  = Z A.  (k)  x (k-a.)  + Z B.  (k)  u (k-b.) 


i=0 


1 j=0  3 


(1) 


where  A^  and  B^  are  respectively  time-varying  n x n and  n x p matrices, 
x is  the  n-dimensional  state  vector,  and  u is  the  p-dimensional  control 
vector  with  p < n. 

a.  and  b.  are  delay  indices  with  a = b = 0,  a.  > a.  ,,  and 
i 3 oo*i-  i-l’ 

b.  > b.  .. 

j - j-1 

The  output  uncertain  observation  can  be  described  by 
y(k)  = Y (k)  C (k)  x (k)  (2) 


where  the  matrix  C is  m x c dimensional  and  y(k)  is  a binary  quantity 
taking  on  values  of  0 or  1 with  the  following  mutually  exclusive 
probability 
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Pr.ty(k)  = 1]  = q(k)  (3) 

and  Pr.[y(k)  = 0]  = 1 - q(k)  (4) 

where  l-q(k)  is  the  probability  of  occurrence  of  a false  alarm,  and 
E [y(k)  Y (j)]  = q(k)  q(j) , k^j  (5) 

E [y2  (k)  ] = q(k)  (6) 


In  some  applications, the  state  variable  x may  not  be  directly 
available,  or  the  dimension  of  the  state  model  is  too  large  that 
it  may  be  desirable  to  reduce  the  order  of  the  estimator.  The  general 
form  of  such  an  observer  can  be  described  by  the  relation 
N N 

z(k+l)  = I F.  (k)  z (k-a.)  + I D. (k)  y (k-a.) 

i=0  1 1 i=0  1 1 

M 

+ Z G (k)  u (k-b.)  (7) 

j=0  J J 

and  the  estimate  of  x takes  on  the  expression 

x (k)  = P(k)  z (k)  + H(k)  y (k)  (8) 

where  z is  i-dimensional  and  n > H > n-m.  F.,  D. , G.,  P and  H are 
of  the  order  jLxt,  , £x  m,  Hx  p,  n x^and  n x m,  respectively. 

Remark  1 : If  there  exists  an  f-x  n dimensional  transformation  matrix 
T such  that 

Ey  [ z (k)  ] = Ey  [T  (k)  x (k)]  (9) 

then, the  necessary  condition  for  x (k)  to  be  an  unbiased  estimate 
of  x (k)  can  be  easily  obtained  as 

P(k)  T (k)  + q (k)  H (k)  C (k)  = In 


CIO) 
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the 

(11) 

(12) 

Theorem  1:  Given  a matrix  T(k)  which  satisfies  the  relation  of 
Eq.  (9),  then  for  the  dynamic  system  of  Eqs.  (1)  and  (2),  one  can 
construct  an  £-order  observer  in  the  form  of  Eqs.  (7)  and  (8). 


Remark  2:  Using  the  relation  of  Eqs.  (7)  and  (9),  one  can  obtain 

Luenberger's  necessary  conditions 

N N N 

T(k+1)  Z A (k)  = Z F (k)  T (k-a.)  + Z q(k)  D (k)  C (k-a.) 
i=0  1 i=0  1 i=0  1 

M M 

and  Z G. (k)  = T (k+1)  Z B.  (k) 

j=0  3 j=0  3 


Proof:  Let  F.(k)  = T(k+1)  A.(k)  P (k-a.)  (13) 

D.(k)  = T(k+1)  A.(k)  H (k-a.)  (14) 

where  P and  H must  be  chosen  to  satisfy  the  constraining  condition 
of  Eq.  (10). 


Then,  making  use  of  Eqs.  (1),  (2),  (7),  (8)  and  (12)  leads  to 

N 

E [z  (k+1)  - T (k+1)  x (k+1)]  = Z T(k+1)  A. (k) 

Y i=0  1 

P(k-ai)  [z  (k-ai)  - T (k-a^  x (k-a^] 

This  implies  that  once  the  initial  conditions  are  selected 

such  that 

z (-a^  = T (-a.)  x (-a^,  i * 0,  1,  2,  ... 

then.E^  (z  (k)]  = E^[T  (k)  x (k)]  for  all  k = 0,  1,  2,  ...  and,  hence, 
the  existence  of  an  i-order  observer. 


3.0  OPTIMAL  MINIMAL-ORDER  OBSERVER  FOR  STOCHASTIC  SYSTEMS 


Consider  a class  of  stochastic  discrete  time-delayed  systems 


given  by 


■MMi  m 


* 
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N M 

x(k+l)  = Z A (k)  x (k-a  ) + Z B.  (k)  u (k-b.)  + £(k)  (15) 

i=0  1 j=0  J J 

and  the  uncertain  observation  is  described  by 

y(k)  = y (k)  C (k)  x (k)  + n (k)  (16) 

where  £ and  n represent  independent  zero-mean  white-noise  processes 
with  known  correlation  matrices  Q and  R respectively.  The  binary 
sequence  y (k)  is  identical  to  that  of  Section  2. 


The  initial  states  x (0)  and  x (-a.),  for  i=0,  1,  ....  N 
are  independent  zero-mean  Gaussian  random  vectors  with  a positive 
definite  covariance  matrix 


Ex  Ey  {X  ('ai5  x T (_ak)}  = Vo  (ai*  V for  i = °»  ....  N and  k 

— 0,  1,  • • . , j . 


Theorem  2:  Given  the  stochastic  system  of  Eqs.  (15)  and  (16)  with  the 
observer  constructed  in  the  form  of  Eqs.  (7)  and  (8),  the  error 
variance  equation  is  described  by 

N 

V£  (k+l)  = [I  - q (k+1)  H (k+1)  C (k+1)]  [ Z A.  (k) 

i,j=0  1 

Vx  (k-ai»  k~aj)  Aj  T (k)  + Q (k)J  (17) 


where  the  optimal  value  of  H (k+1)  minimizing  the  error-variance 
cost  function 


J = Ex  Ey  { x T (k+1)  M ^k+1)  x (k+1i} 


(18) 


is  given  by 

N 

[ ( Z 

.j 


Z A (k)  (k-a  k-a.)  A. 
i,j=0  1 x 1 V J 


00 


+ Q (k)>  C (k+1)  q (k+1)]  [R  (k+1)  + q2  (k+1)  C (k+1) 


H(k+1) 
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{ Z A (k)  (k-a  k-a  ) A T (k)  + Q(k)}  C T (k+1) 
i,j=0  1 J J 

+ { q (k+1)  - q2  (k+1)  } C (k+1)  Vx  (k+1)  C T (k+1)]'1  (19) 

where 

N . - T T 

V (k+1)  = (k+1)  + Z A. (k)  x (k-a.)  x ‘(k-a.)  A.  1 (k) . 

i,  j=0  1 1 3 3 

Proof:  Since  x (k+1)  = x (k+1)  - x (k+1),  utilizing  the  expression 
of  Eqs.  (7)  and  (8)  yields 

N 

x(k+l)  = [I  - H (k+1)  q (k+1)  C (k+1)]  [ Z A.  (k) 

i=0  1 

x (k-a.)  + £ (k)]  + H (k+1)  [q  (k+1)  - Y (k+1)] 

C (k+1)  X (k+1)  + H (k+1)  n (k+1)  (20) 

The  error-variance  equation  is  now  obtainable  by  taking 
the  expectation  of  Eq.  (20)  multiplied  by  its  own  transpose.  The 
optimal  value  of  H (k+1)  can  be  derived  by  setting  the  gradient  of 
Eq.  (18)  equal  to  the  null  matrix  [10].  The  result  is  simply  that  of 
Eq.  (19),  which  in  turn  reduces  the  error -variance  equation  to  Eq.  (17). 

Remark  3:  Substitution  of  Eqs.  (7),  (12)  - (14)  into  Eq.  (8)  also  yields 

N 

x(k+l)  = [I  - q(k+l)  H(k+1)  C (k+1)]  [ Z A.  (k)  x (k-a.) 

" M i=0  1 

+ Z B.  (k)  u (k-b.)]  + H (k+1)  y (k+1)  (21) 

j=0  3 3 

Notice  that  the  estimator  of  Eq.  (21)  together  with  Eqs.  (17) 
and  (19)  constitute  Nahi's  estimator  [7]  for  a linear  discrete  system 
with  uncertain  observations.  In  the  special  case  that  q(k)  = 1,  the 
above  algorithms  reduce  to  Kalman's  optimal  filter.  It  may  also  be 
noted  that  the  optimum  value  of  H(k+1)  and  the  error-variance  equation 
of  Eq.  (17)  are  entirely  independent  of  the  control  vector  u. 
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On  the  other  hand,  the  output  vector  z(k+l)  of  Eq.  (9)  is  an 
estimate  of  T(k+1)  x (k+1)  with  an  observer  "ror  e(k+l)  that  can  be 
easily  found  to  satisfy  the  equation 
N 

e(k+l)  = E [F.  (k)  { z(k-a.)  - T(k-a.)  x (k-a.)  } + D.(k)  n (k-a^ 
i=0 

+ (k)  { Yfk-ap  - q(k-ai)  } C (k-a..)  x (k-ai)] 

- T (k+1)  £ (k).  (22) 

and  the  error-variance  equation  of  e(k+l)  is  given  as 

N 

V£(k+1)  = T (k+1)  [ E A.(k)  P (k-a  ) V (k-a.,  k-a  ) 

i,j=0  J 

PT(k-a . ) A.  T(k)  + A. (k)  H (k-a.)  R (k-a.,  k-a.) 

3 3 l l l 3 

HT(k-a.)  A.T(k)  + Q (k)]  TT(k+l) 

N 

+ E T(k+1)  (q(k-a.)  - q2(k-a.)}  A (k)  H (k-a.) 
i=0 

C (k-a.)  Vx  (k-a.)  C T(k-a.)  HT(k-a.)  A.T(k)  TT(k+l)  (23) 

Now  the  design  of  the  observer  is  to  be  concentrated  on  the 
selection  of  the  matrices  P and  T,  satisfying  the  constraint  of  Eq.  (10). 

In  this  section,  the  observation  equation  of  Eq.  (16)  is 
assumed  to  contain  mj  dimensional  noise-free  components,  namely 

ytkl  - Ytk)  [c‘S]  * M ♦[;,(|0]  M 

where  Cjand  C2  are  respectively,  mj  x n and  m-mj  x n dimensional, 
with  noise  correlation  matrix 

R fkl  = 0,1,1  mi  0 ‘ 

ll  } L 0 Rl  (k)J 

where  R) (k)  > 0.  For  this  case,  it  can  be  easily  shown  that  the 
dimension  of  the  optimal  minimal -order  observer  is  n-m^  [3]. 
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Now,  partition  the  matrix  C(k),  H(k)  and  the  estimator  x(k) 


in  the  forms  of 


C(k) 

H(kj 


Ci i (k) 

C2l(k) 

Hi  (k)l 
_H2  (k)J 


Ci  2 00 
C22(k) 


(25) 


Using  the  compatibility  theorem  of  Tse  and  Athans  [3]  and 
following  the  procedure  similar  to  that  of  Yoshikawa  and  Kabayaski  [8], 
one  obtains 


q(k)  Hj(k)  = Cn'1  (k)  { [I  1 0]  - q(k)  Ci2(k)  H2  (k) } (26) 

and  after  some  algebraic  manipulations,  one  can  show  that  the  pair 


P(k)  = 
and 

T(k)  = 


r-Cn-1(k)Ci2(k) 

L ^n-mj 


jj-q(k)  H2(k)  Cj  (k)  | In  mi  - q(k)  H2 (k)  C2(k)J 


(27) 


(28) 


specifies  one  optimal  minimal-order  observer,  and  Hj (k)  is  uniquely 
specified  in  terms  of  H2(k),  if  q(k)  t 0. 


4.0  GENERAL  REDUCED-ORDER  OBSERVERS 


This  section  considers  a general  class  of  measurement  models 


y(k)  = Y(k) 


(29) 


where  xi  is  mj  dimensional  and  x2  is  n-mi  dimensional,  since  in  many 
physical  systems  some  of  the  measurement  outputs  are  simply  the 
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corrupted  signals  of  the  states  involved.  And  it  is  further  assumed 
that  mi  components  of  the  measurement  vector  are  relatively  noise-free 
compared  to  the  remaining  outputs.  The  problem  here  is  to  construct 
a reduced-order  observer  of  n-mi  dimensional  using  these  mi  measurement 
outputs. 

It  could  easily  be  shown  that  in  the  case  of  y(k)  = 1,  N = 0, 
Ci  * 0,  Eqs.  (15)  and  (29)  constitute  the  canonical  class  of  systems 
considered  by  Leondes  and  Novak  [6].  When  y(k)  = 1,  N = 0,  mi  = m, 
the  solution  of  the  problem  is  reduced  to  that  of  Leondes  and  Novak  [5] 
If  one  were  to  set  y(k)  = 1,  mi  = 0,  N = 0,  the  problem  becomes  the 
special  case  of  Kalman  filteT  [9],  and  when  mi  - 0,  N = 0,  the  solution 
is  equivalent  to  that  of  Nahi  (7). 

Therefore,  the  modeling  of  Eqs.  (15)  and  (29)  does  not  cover 
only  the  reduced-order  observer  theory,  when  all  measurement  signals 
obtained  are  assumed  to  contain  the  state  variables.  It  also  presents 
a rather  broad  spectrum  of  reduced-order  observers,  when  the  detected 
signals  are  associated  with  the  probability  of  false  alarm. 


Based  upon  the  experience  of  Leondes  and  Novak  [5],  [6]  in 
constructing  a classical  form  of  reduced-order  observers  satisfying  the 
constraining  relation  of  Eq.  (10),  one  can  set 


(31) 


From  Eq.  (10),  it  can  be  shown  that 


T(k)  = [-q(k)  (Ki  + K2  C!>  | Inmi  - q(k)  K2  C2] 

Application  of  Eq.  (16)  and  the  definition  of  e(k) 
brings  about 


(32) 


t 
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V£(k)  = [q(k)  - q2(k)]  H (k)  C (k)  Vx(k)  C (k)  H T(k) 

• " »>  ■»”  tv  j ra 

Partition  Eq.  (23)  in  the  form  of 

vk+i>  ■ [;*;S  !'$]]  t T <kti> 

then,  the  trace  of  Eq.  (33)  can  be  derived  using  Eqs.  (30)  - (32). 
From  the  relations  of  Eq.  (7),  (22),  and  the  definitions  of  e(k+l) 
and  W(k) , one  can  show  that 

V (k+1)  = W(k)  + l A.  (k)  x (k-a.)  x T (k-a  ) A.  T(k) 
i,j=0  1 1 3 3 


Finally,  Kj  and  K2  are  optimized  using  the  matrix  minimum 
principle;  the  results  are: 

[Kn  (k+1)  K2 (k+1))  = |q(k+l)Wl2(k)  |q(k+l){W12(k)  C1T(k+l)+W22(k)C2T(k+l)}j 

’Rn(k+1)  + q2(k+l)Wn(k)  + (q(k+l)  - q2(k+l)}  V^^k+l) 

Rl2(k+1)  + q2 (k+1)  { Ci(k+l)Wn(k)  + C2(k+l)W12(k)  } 

+ { q(k+l)  - q2(k+l)  ) { C2(k+l)VXi2(k+l)+Ci (k+1)  Vxii(k+1)  } (35) 

Rl2 (k+1)  + q2 (k+1)  { WU  (k)  CiT  + Wi2  (k)  C2  T(k+1)  } 

+ { q(k+l)  - q2 (k+1)  } { (k+1)  C2  T(k+1)  + V (k+1)  Ci  T(k+1)  } 

R22(k+1)  + q2 (k+1)  Cx (k+1)  { Wn(k)  Ci  T (k+1)  + W12  (k) 

C2  T (k+)  ) + q2 (k+1)  C2  (k+1)  W22(k)  C2  T (k+1) 

+ { q(k+l)  - q2 (k+1)  } { C2(k+1)  V (k+1)  C2  T(k+1) 

x22 

+ Cl (k+1)  V (k+1)  Cl  T(k+1)  + 2 Cl (k+1)  V (k+1)  C2  T(k+1)  } 

X11  x12  J 


4- 
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where  the  noise  covariance  matrix  of  (vj,  V2)  and  the  covariance  of 
x(k+l)  have  been  respectively  partitioned  as 

[Rll(k+1)  R12(k+ljl 

R2l(k+1)  R22(k+l)J 

and 


Vx(k+1) 


rv 

= xn 

V 

Y01 


(k+1)  V (k+1) 

x12 


(k+1)  V (k+1) 
x21  x22 


f 

I. 


In  order  to  guarantee  a unique  solution  for  Kj  and  K2,  it  is 
sufficient  that  Qi 1 and  R22  be  positive  definite. 


The  observer  can  be  initialized  by  defining 
N 

W(0)  = l A (0)  V (a  a ) A.  (0)  + Qn 
i , j=0  3 J 

then  Kj,  K2,  V£,  W and  can  sequentially  be  evaluated  using  Eqs.  (23), 


(33)  - (3S) . 


5.0  ILLUSTRATIVE  EXAMPLE  OF  A RADAR -TRAC KING  PROBLEM 


In  order  to  illustrate  the  performance  of  the  proposed 
reduced-order  observer,  the  radar-tracking  system  of  Leondes  and 
Novak  [6]  is  considered: 


-X(k+ir 

_ T « 

1 T ~2 

’x(k)' 

'0" 

x(k+l) 

= 

0 1 T 

x(k) 

+ 

0 

_x  (k+1 )_ 

.0  0 1. 

.x(k). 

,1_ 

The  measurement  model  is  given  by 
[l  0 0]  fx(k)1 

yW-rOO  [0  j 0J 

where  y(k)  is  a binary  sequence  of  0 or  1 with  probability 
Prob.  [y (k)  = 1]  = p. 


I 


~ ! 
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« 

In  the  numerical  simulation,  three  different  values  of  p were 
selected.  When  p = 1,  the  Kalman  filter  estimates  of  a 30-run  average 
and  the  estimates  of  the  reduced-order  observer  are  compared  to  the 
actual  state  values.  When  p = 0.8  and  0.6,  the  estimates  of  the  reduced- 
order  observers  are  compared  to  those  of  Nahi  estimates.  From  the 
results  of  Figs.  1-5,  it  can  be  seen  that  estimates  obtained  from  the 
reduced-order  observers  are  slightly  less  accurate  than  those  of  Kalman 
or  Nahi  estimates.  As  expected,  the  lowest  error  variances  occur  when 
there  is  a zero  probability  of  false  alarm. 

6.0  CONCLUSIONS 

This  paper  has  extended  the  Luenberger's  observer  to  deter- 
ministic  and  stochastic  discrete  time-delayed  systems,  when  the  detector 
decision  processess  are  associated  with  a non  zero  probability  of  false 
alarm.  The  derivation  makes  use  of  the  matrix  minimum  principle  to 
minimize  the  error-variance  cost  function.  It  can  be  easily  shown 
that  the  reduced-order  observers  of  Leondes  and  Novak  [5],  [6],  Kalman 
[9]  as  well  as  Nahi  [7]  estimates  may  be  considered  as  special  cases 
of  the  proposed  reduced-order  observer  solution.  The  result  is  illus- 
trated by  a radar-tracking  problem. 
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APPENDIX 


A. 


Nahi's  filter  (mj  = 0) 


In  this  case,  the  observer  state  vector  z is  n-dimensional , 


Ki  = 0,  K2  = K,  H = K,  P = I 


and  T(k)  = I - q(k)  K(k)  C(k). 


Substitution  of  Eqs.  (7)  and  (10)  into  (8)  yields 
N 


x (k+1)  = Z A. (k)  x (k-a.)  + K(k+1)  { y(k+l) 
i=0  1 1 

N 


- q(k+l)  C (k+1)  Z A. (k)  x (k-a.)  } 
i=0  1 1 


and  from  Eq.  (34) 


V£  (k+1)  = [q(k+l)  - q2 (k+1) ] K (k+1)  C (k+1) 


Vx(k+1)  C T (k+1)  K T (k+1)  + Ve(k+1) 


+ K (k+1)  R (k+1)  K (k+1) 


One  also  has  the  relation 


N 


W(k)  = Z A. (k)  [V  (k-a., k-a.)  + K(k-a.) 

i , j=o  1 e 1 y 1 

T 


R(k-a^.k-aj)  K (k-a^)  + ( q(k-ap  - q2(k-ap  } 


K(k-a.)  C (k-a.)  Vx(k-a.)  C T (k-a.)  K T (k-a.) 


6k-a.,k-a.]  A.  T (k)  + Q(k) 
N 


Z A (k)  (k-a., k-a.)  A.  (k)  + Q(k) 

i , j=0  1 x 1 J J 


= (k+1 | k) 


(36) 


(37) 


(38) 


j 
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and  Eq.  (35)  becomes 


K(k+1)  = q(k.+l)  V^(k+l|k)  C 1 (k+1)  R(k+lJ 

+ q2(k+l)  C (k+I)  VA.(k+l|k)  C f (k+1) 

X -1 

{ (q(k+l)  - q2 (k+1)  } C (k+1)  Vx(k+1)  C T (k+1)  } (39) 

Finally,  application  of  Eqs.  (33),  (37)  and  (39)  leads  to 


V£(k+1)  = [I  - q(k+l)  K (k+1)  C (k+1)]  V£(k+l|k)  (40) 

It  is  now  obvious  that  Eqs.  (36)  to  (40)  constitute  the 
Nahi  type  filter  under  uncertain  observations,  and  in  the  special 
case  of  y (k)  = p(k)  = 1,  they  are  reduced  to  Kalman  filter. 

B.  The  equivalence  of  the  results  presented  in  this  paper  and 
those  of  Leondes  and  Novak  [5],  [6]  can  easily  be  established  using 
Eqs.  (23),  (30)  - (33)  and  (35). 
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and  Athans  and  Leondes  and  Novak's  reduced-order  observers.  The  design  and  Athans  and  Leondes  and  Novak's  reduced-order  observers.  The  design 

of  a reduced-roder  observer  is  illustrated  using  a radar- tracking  of  a reduced-roder  observer  is  illustrated  using  a radar-tracking 

problem.  (0)  problem.  (U) 
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